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ABSTRACT 

This report describes a method for analytically determining 
the amount of twist and lean of a tall tower of equilateral 
triangular cross section. This method is also applicable to tall 
structures of other shapes and cross-sectional areas. 

A BASIC computer program that determines the angle of 
twist and amount of lean is provided in Appendix A. 
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SECTION I 
INTRODUCTION 


1.1 PURPOSE 

This report presents an analysis of recorded data for the 
translational displacement of the center of mass and the rotational 
motion about the center of mass for the horizontal section of a 
tall tower. 

1.2 SCOPE 

This report contains illustrations showing a tall tower, 
diagrams of displacement measurement geometry and mechanics, and 
an analysis of translational and rotational displacement data. 

1 . 3 BACKGROUND 

A 150-meter tower at the Kennedy Space Center (KSC) is 
presently used as a structure for gathering wind data (Figure 1-1) . 

In the course of this operation, the question arose as to the 
amount of tower twist and lean. Preliminary measurements established 
that there was a slight twist in the tower. Simultaneous measure- 
ments on the three vertices of the tower at given cross-sectional 
levels (Figure 1-2) were taken to determine the displacement of these 

vertices. By the method described in this report, these measure- 
ments were reduced to show the amount of twist and lean of the 
tower. 

1.4 ACKNOWLEDGEMENTS 

Acknowledgement for their contribution to this report is made 
to the personnel of LV-GDC-1 Gyro and Stabilizer Systems Branch 
for the simultaneous theodolite measurements. 
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Figure 1-1. 


Overall View 
Ground Winds 


of 150-Meter 
Tower 
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SECflON II 
ANALYSIS 


2 . 1 GENERAL 

The displacement of any rigid body can be described as 
the translational displacement of the center of mass and the 
rotational motion about the center of mass. 

In this analysis, each horizontal section of the tower at 
the respective measurement level was treated as a symmetrical 
rigid body free to be translated or rotated only in a horizontal 
plane. Figure 2-1 is a schematic representation of such a hori- 
zontal section, and Figures 2-2 and 2-3 show the measuring 
geometry and mechanics. 



Figure 2-1. Schematic Representation of a 

Horizontal Section of the Tower 

2.2 MEASUREMENT GEOMETRY AND MECHANICS 

A theodolite was positioned at a range, r, from the tower 
at the intersection of the arrows (point E, Figures 2-2 and 2-3). 
Observations were made at the various levels of the tower, and 
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measurements were taken from the inside edge of the vertex strut 
at each level (point B, Figures 2-2 and 2-3) using the 67-foot 
(20 .4-meter) level as the point of zero deflection. The theodo-' 
lite was elevated and the angular displacement of the vertex strut 
at successive levels was recorded. 

First, a relationship had to be established between the mea- 
sured displacement at the edge of the tower and the associated 
relative displacement at the center of the side toward the theodo- 
lite (point E, Figures 2-2 and 2-3), due to a twist in the tower. 
The difference in measured displacements due to any translational 
motion of the tower is negligible as long as r is very much greater 
than i. For proof, consider Figure 2-2, where the following nota- 
tion applies : 


C' center of tower 

2i length of side of tower (=2BF) 

Ai measured displacement at level in question (=BD) 
r ground distance of theodolite from tower (=EF) 

h relative displacement at center of tower side (=FG) 


If C moves to C' in a pure translation, then point F will 
move to G and point B will move to D. Since triangle BCF equals 
triangle DC’G, CC’=FG, and h=Ai. 4 


now, tan /FEG = — 

r 

ZFEG - arc tan — 

r 

Also, tan /FEB - - 

r 

tan /FED = 1 tAl = A.+. h 

r r 

Therefore, ,/FED - /FEB = /BED = arc tan - +h - 

r 


Now angle BED is the observed angle, while the 
angle is FEG . The difference between these angles is 

A = arc tan - ~ - arc tan ^ - arc tan — 
r r r 


desired 


Expanding the arc tangent functions, we obtain: 

A = - I il *h)3 + 1 (i +h) s 

r 3 r 3 5 r 5 

- i + I ii 1 if. 

r 3 r 3 5 r 5 

_ h + 1 h3_ 1 h 3 

r 3 r 3 5 r 5 
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Figure 2-2. Translational Displacement Measurement 
Geometry and Mechanics 
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Figure 2-3. Rotational Displacement Measurement 
Geometry and Mechanics 
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All terms involving h and i alone vanish, leaving: 
1 3i 2 h+3ih 2 


A = - 


+ I 5i 4 h+10i 3 h 2 +10i 2 h 3 +5ih 4 
5 r 5 


The resulting error in the measurement of h would be: 


A- r= - 

i 2 h+ih 2 

L i 4 h+2i 3 h 2 +2i 2 h 3 +ih 4 

r 2 

r 4 . 

Let h 

= 1, i = 

4, r = 400, then: 


20 

422 

A • r= - 

+ 

— — — = -1. 25xl0" 4 


16xl0 4 

256xl0 8 

A« r=- 

1.25x10- 

4£_ 1 

TO*o 


This is the resultant error when h is subsituted for Ai in the 
example . 

In Figure 2-3, we have the following notation: 

center of tower 
length of side of tower 

measured displacement at level in question (=BD) 
relative displacement at center of tower side (=FG) 
derived displacement at center of tower side f=IH) 

(used later in analysis) 
distance from center to vertex (=BC) 
distance from center to side (-CF) 
twist angle ’ 

ground distance of theodolite from tower (=EF) 

The mechanics of Figure 2-3 are explained as follows. As 
tower twists, the line CF moves through an angle $ to position 

LH. The diagonal CB moves through the same angle $ to position CA 

Since the theodolite measures angles, the instrument sees a dis- 
placement from point B to point D, which is the point where the line 
AE intersects the extended line of the tower side. If the theodolite 

had a reference point F at the center of the tower, it would see a 

displacement from point F to point G which is the intersection of HE 
with the tower side. First, a relationship must be found between 
BD, i.e. ai, and t, i.e. s sin $ ( = IH) , which latter quantity is 
V S cr\ in the anal y sis » an< * finally a relationship between t and h 


C 

2i 

Ai 

h 

t 

d 

s 

$ 

r 
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From Figure 2-1 we have the following relationships: 
cos 30° = tan 30° » j 

d = — = s = i tan 30° = — d = 2s 

cos 30° /3 /3 

From Figure 2-3, 

t = s sin $ 

From the law of cosines, triangle ABC, 
a 2 = 2d 2 - 2d 2 cos $ = 2d 2 (1-cos $) 

a = /2d* /1-cos $ = /2d / 2(1 ~ cos 
a = 2d sm ^ = 4s sm y 


Also , 

2 /CBA + 4> = 180° 
/CBA = 90° - | = a 
In triangle ABD, 


/ABD + /CBA + 30° = 180° 

/ABD = 180 - 30°-/CBA = 180° - 30° - 90° + | 


/ABD = 60° + | = & 


In triangle ACE, from the law of tangents, 

tan C^CAD - /ABC) = r+s-d (60° + ») 

2 r+s+d 2 

/CAD -/ABC = 2 arc tanjf^f^ cot ^ + = 2 a 

I- 

/CAD -ZAEC = 2K, where arc tan 


r-s 60° + * 

r+3s cot 2 


arc tan 
= K 


r-s 60' 

r+3s cot 


+ $ 

T~ 
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Again in triangle ACE, 

ZCAD + ZAEC + 60° + $ = i80° 
ZAEC = 120° - 4- -ZCAD 
ZCAD - 120° + * +ZCAD = 2K 
2 ZCAD = 2K + 120° - $ 

ZCAD = K + 60° - | 


Also, 


ZDAB = ZCAD - a = K + 60° - | - 90° + - 

2 2 

ZDAB = K - 30° 

ZADB + g +ZDAB = 180° 

ZADB = 180 - g - K + 30° = 180° - 60° - | - K + 30° 
ZADB = 150° - | - K 

In triangle ADB, by the law of sines, 


sin ZADB 


Ai 

sin ZDAB 


i = a sin (K- 30°) 

sin (150° - 1 = K) 


Ai = 4s sin ^ 


sin (150* 


5 

2 

Now t = s sin 4> = s sin 2 -| 


sin ('K-30°‘) 


- T - K) 


t = 2s sin ~ cos | 
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2t 


4s sin y cos j 


4 s sin 



2t 


<f 

cos 2 


Substituting in the expression for.Ai, 


II 

♦H 

< 

2t 


sin 

(K- 30° ) 


$ 

COS — 
2 


sin 

(150° - | - 

K) 

t_ = 

I COS 

<i> 

sin 

$ 

(150° - J - 

ja 

Ai 

2 

sin 

(K - 30°) 



Figure 2-4 is a plot of versus $ , using the following known 
quantities: 1 

r = 400 feet 
i = 3.833 feet 
s = 2.213 feet 
d = 4.426 feet 

Use will be made of this graph (Figure 2-4) later. (It will be 
noted that for $ = 0, the problem is not defined.) 

In order to determine how much t differs from h, we proceed as 
follows : 

In triangle FEG , tan ^FEG = ^ 

also, tan ^FEG = tan ^IEH = r . (s ! s c05 „■> 

h = t 

r r + Cs-s cos 4>) 


h _ r 

t r + s (1-cos 4>) 


Now , 


sin 2 


$ _ 1-cos $ 
2 2 


therefore , 


h 

t 


r 

r + 2s sin 2 j 
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If r = 400 feet 
s = 2.213 feet 
$ = 10 degrees 

We have 

h _ 400 _ 400 

t 400 + 4.426 T- 0075969) “ 400.03362 

~ = 0.99992 


Therefore in this case, h differs from t by 1 part in 12,000. 


2.3 MEASUREMENT ANALYSIS 

We can now determine the twist angle and the displacement 
of a given level of the tower from the three thodolite measure- 
ments. Let us first define a quantity Ai r , the actual measured 
displacement, which is the sum of the displacements due to trans- 
lation and rotation. Also we define a quantity t' which is the 
sum of the t's due to translation and rotation. In other words, 
each Ai 1 can be considered as being the sum of three components: 
a component due to the x displacement, a component due to the y 
displacement, and a component due to twist. In like manner, each 
t' can be considered as being made up of three components: a 

component due to the x displacement, a component due to the y 
displacement, and a component due to twist. The components due 
to the horizontal displacements are the same, (with negligible 
error) for both Ai’ and t’, while the components due to twist 
are related by the function graphed in Figure 2-4. 

Now consider Figure 2-5 where the total derived displace- 
ment, that is, the t’s for each theodolite are shown. The x 
displacements are positive to the right, y displacements are 
positive upward, and angles are positive clockwise. We may write 
the following two sets of equations: 

t| = 0«Ax + Ay + s sin $ 

t£ = -Ax sin 60° - Ay sin 30° + s sin 4> 

t| = Ax sin 60° - Ay sin 30° + s sin * 

Aij = 0 ‘Ax + Ay + $ ($) ♦ s sin $ 

Ai 2 - -Ax sin 60° - Ay sin 30° + $ ($■) * s sin $ 

A i 3 = Ax sin 60° - Ay sin 30° + $ ($) • s sin $ 

where$ ($) stands for some function of $. 
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Figure 2-5. Tower Displacement Components 


2-11 





TR- 297 
Rev 1 




The second set of equations constitute three simultaneous 
equations in three unknowns, and the solution is readily obtain- 
able by determinants. The unknowns are Ax, Ay, and sin $>, and 
the known terms are AiJ, Ail', ...and Ail. The determinant D of the 
coefficients of the unknowns is: 



0 

1 

S 


0 

1 

1 

D = 

-/T/2 

-0.5 

S 

= S 

-/T/2 

-0.5 

1 


+ /T/2 

-0.5 

S 


+ /3/2 

-0.5 

1 


Therefore, Ax can be written: 



Ai i 1 S 


Ai 1 

1 

1 

A 1 

Ax = B 

Ail -0.5 S 

S 

D 

Ail 

-0.5 

1 

where S ! 

Ai^ -0.5 S 

= *(*) • s 


Ail 

-0.5 

1 

Ax = ITT 

(-0.5 Ail + 0.5 Ail 

+ 0.5 Ai J + 

Ail 

0.5 Ai{ ~ 

Ail) 


Ax = -3 a i 2 + 3 AiJ = Ai^ - Ai£ 

3/T 7 ? 

Ax = & ( Ai 3 - Ail) 

Likewise for Ay: 


0 

Ail 

S 


0 

Ai{ 

1 

- /T/2 

Ail 

s 

s 

D 

-/T/2 

Ail 

1 

+ /T/2 

Ail 

s 


+ /3/2 

Ail 

1 


Ay j ^ Ail ■ /T/2 Ail + /T/2 Ai{ + /T/2 Ai{) 


Ay 


2 Ai{ - A i 2 - A i 3 
3 
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And for $ O) sin 0 , 

0 1 Aij 

* (*) sin <}> = -/T/2 -0.5 Ais 

+/3/2 - 0.5 A i 3 

* (*) sin 4> = 2_ (/T/2 Ai .3 + /I/2 Ai£ = /T/4 Ai{ + /T/4 Ai[ ) 

3s /T 


$ ($) sin $ 


Ai[ + Ai£ + A i 3 


sin $ - f ^ (Ai l + Ai 2 + Ai 3) 
3s 


where f f$) = — — — 

* O) 

+ A i^ "*■ Ai 3 ) 

For the NASA 150-Meter Ground Winds Tower, s = 2.213 feet, and 
therefore: 

Ax = 0.577 (a i 3 - Ai^) 

Ay = 0.333 (2 A i { - Ai£ - Ai^) 

$ - arc'sin [0.01255 f ($) (Ai{ + Ai| + a i 3 3 ] 
where the coefficient 0.01255 is in inches" 1 . 



It will be noted that the correct value of the twist angle, 
must be known before f ($) can be found. However f (<j>) can be 
readily found from Figure 2-4 by successive approximations, assum- 
ing an initial,. arbitrary , value of <t> . An example will make the 
method clear. Some actual tower measurements are shown in Table 
2 - 1 . 
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Table 2-1. Tower Twist and Lean Measurements 


Level 

Altitude (ft) 

Ai{ (in.) 


| 

1 

1 (Ref) 

67.10 

0 

“ — - —i 

0 

0 

2 

133.70 

-1.08 

-0.24 

-0.94 

3 

192.10 

-1.49 

-0.40 

-2.12 

4 

250. 40 

-1.51 

-0.19 

-1.93 

5 

308.70 

-1.95 

+ 0.21 

-4.10 

6 

367.10 

-2.66 

+0.09 

-5.30 

7 

425.40 

-3.30 

+0.66 

-6.22 

8 

483.73 

-3.71 

I 

+ 0.85 

-7.11 


assume a value for $ of 5 degrees at level 8. From 
the graph of Figure 2-4, we obtain a value for f ($) of 1.104. 
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This will give a new value for $ of, 

♦ = arc sin (0.01255 x 1,104) (-3.71 + 0.85 - 7.11) 

$ = arc sin (0.013855) (-9.97) 

4 - arc sin (-0.13813) 

4 = -7°56 , 23" = -7.94° 

Again entering the graph of Figure 2-4 we obtain a conversion 
factor of 1.162, therefore, 

4 - arc sin (0.01255 x 1.162) (-9,97) 

4 = arc sin (0.014583) (-9.97) 

4 = arc sin (-0.14539) 

4 = -8°21 ’ 36" = -8.36° 

Repeating the preceding process, 

4 = arc sin (0.01255 x 1.171) (-9.97) 

4 = arc sin (0.014696) (-9.97) 

4 =* arc sin (-0.14652) 

4 = -8°25'31" = -8.42° 

A final iteration gives the required value, 

4 = arc sin (0.01255 x 1.172) (-9.97) 

$ = arc sin (0.014709) (-9.97) 

$ = arc sin (-0.14665) 

4 = - 8° 25 ’ 5 8 M = -8.43° 

Table 2-2 shows the calculated values of the displacements and 


2-15 



TR-297 
Rev 1 


Table 2-2. Tower Displacement and Twist Angle 
Calculated Values 


Level 

Ax (in.) 

Ay (in.) 

$ (degrees) 

1 

0 

0 

0 

2 

-0.40 

-0.33 

-1.69 

3 

-0.99 

-0.15 

-3.08 

4 

-1.00 

-0.30 

-3.02 

5 

-2.49 

0 

-4.61 

6 

-3.11 

-0.04 

-6.41 

7 

-3.97 

-0.35 

-7.34 

8 

-4.59 

-0.39 

-8.43 


As a check on the calculations, a graphical overlay model 
was constructed of the top level and the ground level. The twist 
angle as measured with a protractor was 8.5 degrees, and the dis- 
placement of the tower center was measured as 4.68 inches as 
against the calculated 4.61 inches. These results are in ex- 
cellent agreement. 

The analysis of the last section is directly applicable to 
structures of circular cross section. If three fiducial marks are 
equispaced around the structure (point F, Figure 2-6, for one mark), 
then the displacement and twist of any cross section may be deter- 
mined from the preceding equations for Ax, Ay, and $. It is 
assumed, in this case, that the theodolite is pointed normally to 
the tangent at the fiducial mark. 
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Figure 2-6. Structure of Circular Cross Section 
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APPENDIX A 
COMPUTER PROGRAM 

A program has been written in BASIC for use with the GE-635 
computer time sharing system. The following notations have been 
used m the program: 


Ail = Ij 

Ai^ = I 2 

Ai 3 - I3 

* “ P (initial value) 

K = FNA (P) 

f(*) = FNC (P) 

Ax « X 

Ay = Y 

$ = P3 (final value) 

LEVEL = L 


CQ T l ie va :[ ue f used in the data steps are the values used in 
the sample calculations. Any values of course could be used. 
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010 REM THIS PROGRAM DETERMINES THE ANGLE OF TWIST AND AMOUNT OF LEAN 
020 REM IN HORIZONTAL X § Y OF A TRIANGULAR TOWER WHERE S IS THE 
030 REM DISTANCE FROM THE CENTER TO ONE SIDE. THREE THEODOLITES 
040 REM AT THE SAME DISTANCE R FROM THE CENTER OF EACH SIDE MEASURE THE 
050 REM HORIZONTAL DISPLACEMENTS II, 12, $ 13 OF THE CORNERS OF THE 
060 REM TOWER RELATIVE TO SOME REFERENCE LEVEL. (SEE KSC TR-297, REV 1} 
070 PRINT SPC(O) ’'LEVEL'’ SPC(9)"X" DISPLACEMENT” SPC(9)”Y” DISPLACEMENT”- 
080 PRINT SPC (9) "TWIST ANGLE" 

090 PRINT SPC (2) "#"SPC (18) "IN . "SPC (20) "IN . "SPC (17) "DEG." 

100 READ R,S 

110 READ L, II, 12, 13 

120 DEF FNA(P) = ATN ( (R-S) *COT ( (60* . 01 74533+P) / 2 ) / (R+3*S) ) 

130 DEF FNB(P) = . 5*COS (P/2) *SIN ( 150* . 01 74533-P/2 - FNA (P) ) 

140 DEF FNCCP1 = FNB (P) /SIN (FNA (P) -30* . 0174533) 

150 REM AT THIS POINT AN ESTIMATE IS MADE AS TO THE APPROXIMATE VALUE OF 
160 REM P, THE TWIST ANGLE. IN THIS CASE P IS ESTIMATED TO BE 5 

170 REM DEGREES, OR 0.87266 RADIAN. 

180 LET P = .087266 

190 LET F = FNC(P) 

200 REM PI = ARC SIN (F* ( I 1+1 2+1 3) / (3*S*12 ) ) 

210 REM NO ARC SIN FUNCTION IN BASIC. 

220 REM ARC STN(G) = ATN (G/SQR(l-G+2) ) 

230 LET G=F*(I1+I2+I3)/(3*S*12) 

240 LET P1=ATN (G/SQR(l-G+2) ) 

250 LET P2=ABS (PI) 

260 IF P=P2 THEN 290 
270 LET P=P2 
280 GO TO 190 

290 LET X=(1/SQR(3))*(I3-I2) 

300 LET Y= (1/3)* (2*11- 12- 13) 

310 LET P3=P1*57 .2958 

320 PRINT USING 330,L,X,Y,P3 

** ###.## ###.## ### 


340 GO TO 110 
350 DATA 400,2.213 

^7n 08 »' ■ 24 ’~ -94,3,-1.49,- .40,-2.12,4,-1.51,- 

lln 21, -4. 10, 6, -2. 66 , + .09, -5. 30 

389 5H A 7 >~ 3 * 39 > + *66,-6.22,8,-3.71,+.85,-7.11 

390 END 
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